12 years of X-ray variability in M31 globular clusters, including 8 
black hole candidates, as seen by Chandra 
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ABSTRACT 

We examined 134 Chandra observations of the population of X-ray sources associated with 
globular clusters (GCs) in the central region of M31. These are expected to be X-ray binary 
systems (XBs) , consisting of a neutron star or black hole accreting material from a close compan- 
ion. We created long-term lightcurves for these sources, correcting for background, interstellar 
absorption and instrumental effects. We tested for variability by examining the goodness of fit for 
the best fit constant intensity. We also created structure functions (SFs) for every object in our 
sample, the first time this technique has been applied to XBs. We found significant variability 
in 28 out of 34 GCs and GC candidates; the other 6 sources had 0.3-10 kcV luminosities fainter 
than ^2x10"^^ erg s"^, limiting our ability to detect similar variability. The SFs of XBs with 
0.3-10 keV luminosities ~2-50xl0'^^ erg s~^ generally showed considerably more variability than 
the published ensemble SF of AGN. Our brightest XBs were mostly consistent with the AGN 
SF; however, their 2-10 keV fluxes could be matched by <1 AGN per square degree. These 
encouraging results suggest that examining the long term lightcurves of other X-ray sources in 
the field may provide an important distinction between X-ray binaries and background galaxies, 
as the X-ray emission spectra from these two classes of X-ray sources are similar. Additionally, 
we identify 3 new black hole candidates (BHCs) using additional XMM-Newton data, bringing 
the total number of M31 GC BHCs to 9, with 8 covered in this survey. 



Subject headings: x-rays: general 
individual — black hole physics 

Introduction 



x-rays: binaries — globular clusters: general — globular clusters: 



When studying the X-ray populations of exter- 
nal galaxies, globular cluster X-ray sources are of 
particular interest, because they are most prob- 
ably low mass X-ray binaries (LMXBs), rather 
than e.g. background galaxies, supernova rem- 
nants, or high mass X-ray binaries (HMXBs). 
Globular clusters contain ^10% of the Galac- 
tic LMXB population (14 out of --150), but 
are thought to contain only 0.1% of the mass, 
hence the LMXB density is ^100 times higher in 
GCs than elsewhere (see the comprehensive re- 



view by IVerbunt fc Lewinl l2006i and references 
within); indeed, it has been speculated that all 
LMXBs were born in globular clusters (see e.g. 
iD'Antona et al.ll2007L and references within). 

M31 is th e nearest neighboring spir al galaxy, 
at ~780 kpc (jStanek fc Garnavichlll998l ). and is a 
favored target of many X-ray telescopes. Hence, 
its GC X-ray sources have been particularly well 
studied. X-ray su rveys of M31 include thos e 
made wi th Einstein ( Trinchieri fc Fabbiandll991 ) , 
ROSAT (ISupper et al.lll997Ll200lh XMM-Newton 



(JTrudolvubov fc Priedhorskvl2004IShaw Greening et al 



2009|: | Stiele et al.ll201l|). and Chandra (iKong et al. 



200a bi Stefano et al.l2"ooilWilliams et al.l20ol )^ 



Di Stefano et al.l (|2002l ) conducted a GC survey 
of M31 with Chandra, covering '^2560 arcmin^. 
They found 30 GC X-ray sources, 15 of which were 
new discoveries. They found that the M31 GC X- 
ray population was systematicahy more luminous 
than the Milky Way (MW) population; 10 M31 
GC X-ray sources in their survey exceeded 10"^^ 
erg s~^, vs. 1 of the 11 MW GC X-ray sources 
surveyed by ROSAT (fv^rbunt et al.lll995[) . 

Recently, iPeacock et al.l (|2010l ) conducted a 
surv ey of X-ray source s in the 2XMMi cata- 
log (jWatson et al.l l2009l) that included all pub- 
licly available XMM- Newton observations of M31, 
searching for associations with GCs; the observa- 
tions cover ~80% of the GCs in M31. They found 
41 GC associations in the 2XMMi catalogue , 
and 4 transient X-ray sources observed by either 
ROSAT or Chandra, for a total of 45 GC LMXBs. 

We have been monitoring the central region of 
M31 with Chandra for the last ~12 years, aver- 
aging '-~^1 observation per month, looking for X- 
ray transients; we exclude periods where M31 is 
too close to the Sun to be observed. These ob- 
servations have no fixed roll angle, and the com- 
bined field may be approximated by a circle with a 
^20' radius. The earlier Chan dra surveys included 
long-term variability studies: iKong et al.l (J20021 ) 
found 204 X-ray sources in the central 17' x 17' re- 
gion, detecting variability in '^50% of the sources 
over 8 ACIS-I ob servations between 199 9, October 
and 2001, June; IWilliams et al.l (|2004l ) examined 
2.5 years of HRC data, and found at least 25% of 
166 sources to be significantly variable. 

In this work we present a variability survey 
of all X-ray sources in our Chandra observations 
that are associated with objects in the Revised 
Bologn a Catalog of M31 globular clusters (RBC 
V.4. Galleti et al.ll2004l [2006. 2007. 2009). This 
inc ludes 30 out of 45 GC X-ray sources identified 
bv IPeacock et al.l (|2010l ): the other 7 RBC asso- 
ciations in our survey include 4 GC candidates, 2 
background galaxies (AGN), and 1 star. We com- 
pare the variability of these X-ray sources with 
th e ensemble vari a bility for active galaxies found 
bv lVagnetti et al.l (|201ll ). 

It has long been known that AGN may vary 
by a factor of a 2-3 on time-scales of months 



to years, with the amplitude of variation in- 
versely proportional to the luminosity (see e.g. 
Marshall et"ai]ll98ll : iNandra et aI1ll997l and ref- 
erences within). However, it has been known for 
exceptional AGN to flare up by an order of mag- 
nitude fe.g. lTananbaum et al.lll978l ). 

Recently IVagnetti et al.l ( 20111 ) studied the en- 
semble variability over time-scales of hours to 
years of AGN in the se rendipitous source cata- 
logs fr om XMM-Newton (IWatson et al. l l2009[ ) and 
Swift O Puccetti et al.ll2011 ): their sample covered 
redshifts ~0.2-4.5, and 0.5-4.5 keV luminosities 
^1043-10^6 erg s-i. They included 412 AGN 
From the XMM-Newton catalog, and 27 AGN 
from the Swift catalog; all of these AGN were sam- 
pled at least twice. They used a structure function 
(SF) to estimate the mean intensity deviation for 
data separated by time r: 

SF (r) s y|(|log/x(i + T)-log/x(t)|)'-a2, 

(1) 
where CTh is the photon noise and /x is the X-ray 
flux. They grouped the SF into logarithmic bins 
with width 0.5; each bin in the range log(T) = 
1.0-3.0 contained more than 100 measurements. 



Vagnetti et al.l (|201l[ ) found that SF increased 



from ~0.1 to ~0.2 as r ranged over 0.1-1000 days 
in the rest frame of the AGN. SF -- 0.15 for t = 
30 days, suggesting that typical AGN would vary 
by '^25-75% between our observations. They also 
investigated the well-observed anti-correlation be- 
tween intensity variability (/var) and luminosity, 
expressed in the form /var oc L^''; k ^-^0.3 in the 
literature for time-scales of days to tens of days. 
They measured k for AGN grouped logarithmi- 
cally over log Lx = 43.5-45.5 for two values of r: 
1 day and 100 days. They find k = 0.42±0.03 for 
T = 1 day, and k = 0.21±0.07 for r = 100 days. 



Vagnetti et al.l (|201l[ ) measured the average 



values of a'^ for the XMM-Newton and Swift sam- 
ples to be 0.031 and 0.163 respectively. When 
they used these values to calculate the SFs for 
each sample, the two SFs were consistent; if SF(t) 
oc r^ then b = O.lOiO.Ol for the XMM-Newton 
SF, and b = 0.07±0.04 for the Swift SF. 

We have created structure functions for each of 
our targets from their 0.5-4.5 lightcurves. The 
goal was to ascertain whether we may distin- 
guish between XBs and AGN from their variabil- 



ity alone; this is important because XB and AGN 
often exhibit similar spectra. 

Our analysis of the full data set has allowed 
us to identify seven black hole candidates associ- 
ated with globul ar clusters; four were previously 
published (B arnard et alj 1201 lal) . and three are 
new to this work. We identified one more GC 
black hole candidate iii this region in earlier work 
( Barnard fc Kolbl 2009), a nd one away from the 
bulge ([Barnard et al.ll2008l) . They were all identi- 
fied from their X-ray properties, as they exhibited 
characteristic low state emission spectra at lumi- 
nosities that exceeded the limit for a neutron star 
accretor; oth er neutron star emissio n models were 
rejected (see IBarnard et al.l l2011ai for the most 
complete justification of this method). 

We describe the data analysis in the next sec- 
tion, followed by a presentation of our results, then 
a discussion of our findings and conclusions. 

2. Observations and data reduction 

We have analyzed 89 ACIS observations and 45 
HRC observations, in order to discern the variabil- 
ity of X-ray sources associated with GCs in this 
region. The merged ACIS image was registered to 
the Field 5 B band image from the Local Group 



Galaxy Survey observations of M31 (|Massev et al. 



I2OO6I ) with IRAF software, using 27 X-ray bright 
GCs. After refining the GC positions in the Chan- 
dra and LGS B band images with imcentroid, we 
found the sky co-ordinates of the new positions 
with XY2sky, then registered the X-ray positions 
to the LGS B band image using CCMAP. The mean 
r.m.s. uncertainty in position due to registration 
was 0.11" in R.A. and 0.09" in Dec. 

We obtained 0.3-7.0 keV lightcurves and spec- 
tra from circular source and background regions 
for each source. The background region was the 
same size as the source region, and at a similar off- 
axis angle. The extraction radius varied between 
sources, because larger ofi^-axis angles resulted in 
larger point spread functions. 

For ACIS observations, we obtained corre- 
sponding response matrices and ancillary response 
files. Source spectra with >200 net source counts 
were freely fitted with absorbed power law mod- 
els. For GCs with multiple freely fitted spectra, we 
obtained the best fit constant values for absorp- 
tion (A^h) and photon index (F). For observations 



with <200 net source photons, we fixed Nn and 
F to the best fit values; if a particular GC had 
no freely fitted spectra, then we assumed A^h = 
7x10^° H atom cm"^, and F = 1.7. We obtained 
luminosities for faint sources by determining the 
unabsorbed 0.3-10 keV luminosity equivalent to 
1 count s~^; multiplying the source intensity by 
this conversion factor gave the source luminosity, 
after correcting for the exposure, vignetting and 
background. 

For HRC observations, we included only PI 
channels 48-293, thereby reducing the instrumen- 
tal background. We used the WebPIMMS tool 
to find the unabsorbed luminosity equivalent to 1 
count s~^, assuming the same emission model as 
for the ACIS observations with <200 photons. We 
created a 1 keV exposure map for each observa- 
tion, and compared the exposure within the source 
region with that of an identical on-axis region, in 
order to estimate the necessary exposure correc- 
tion. We multiplied the background subtracted, 
corrected source intensity by the conversion to get 
the 0.3-10 keV luminosity. 

We created long term 0.3-10 keV lightcurves 
for each source, using the luminosities obtained 
from each observation as described above. We 
only included observations with net source counts 
> after background subtraction. We fitted each 
long term lightcurve with a line of constant inten- 
sity, in order to ascertain the source variability. 
If the null hypothesis probability for the best fit 
to a lightcurve was <0.003, then the variability 
of that lightcurve had a significance >3cr, and we 
classified the target as variable; the null hypothe- 
sis probability determines the likelihood that ob- 
served deviations from the model are due to ran- 
dom fluctuations in the data. 

Luminosity uncertainties for freely fitted spec- 
tra were estimated with XSPEC by calculating a 
range of fiuxes obtained by varying the emission 
parameters; the uncertainties for the faint spectra 
are derived directly from intensity uncertainties 
using the relevant fixed model. 

We derived SFs from the 0.5-4.5 keV fluxes of 
each observation of every target by assuming a 
power law spectrum with the same photon index 
as for the HRC and faint ACIS observations; an 
M31 X-ray source with a 0.3-10 keV unabsorbed 
luminosity of 1.0x10^^ erg s~^ and F — 1.7 has a 



Vagnetti et alj (|201l[ ) calculated the noise com- 
ponent from 



a2=2((51og/x)')^2(loge) 




(2) 

assimiing that Sfx/fx = (l/A^phot)°'^, and iVphot 
is the number of photons. Our lightcurves are 
background subtracted, and ARF-corrected; fur- 
thermore, uncertainties in the luminosities of 
bright sources include uncertainties in the spectral 
parameters. As a result, our uncertainties are not 
simply due to photon counting noise. Hence in 
our case. 
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The merged observations cover an approxi- 
mately circular region with 20' radius. We de- 
tected 430 X-ray sources in this region, and found 
428 globular clusters from the RBC. Therefore 
we expect chance coincidences of X-ray sources 
within 1" of the GC centers for 0.12 out of 428 
GCs. 

Shifting the declination of each source by ±5" 
resulted in no coincidences within 1" of a RBC 
object; shifting the RA by —5" also resulted in no 
coincidences; however, shifting the RA by -1-5" re- 
sulted in two matches within 1" of a RBC object. 
Expanding the search radius from 1" to 2" pro- 
vided one additional match each for offsets ±5" in 
declination. Hence we expect 0-2 false matches. 

We found 37 X-ray sources to be coincident 
with M31 GCs, as classified by the RBC. X-ray po- 
sitional uncertainties include uncertainties in the 
centroid on the X-ray image, uncertainties in the 
registration (0.11" for RA, and 0.09" for Dec), and 
0.25" uncertainty in the LGS position. 

In Table [1] we present the positions and iden- 
tifications of each X-ray source that we associate 
with a member of the RBC V4. We provide the 
latest classific ations of the optical counterparts, 
gleaned from I Caldwell et al.l (|20Q9l ): our sample 
includes 1 star, and 2 galaxies, as well as 30 con- 
firmed old GCs and 4 GC candidates. We also 
show the distance between the X-ray source and 



its optical counterpart, and the uncertainties in 
RA and Dec; X-ray sources with no good centroid 
solution are indicated by "...". Eight of our tar- 
gets are far off- axis (~ 11-20') and also X-ray faint. 
We estimated the positional uncertainties for these 
targets by binning the merged image by a factor 
9 in a; and y, then running imcentroid on this 
binned image. 

Tabled] details the X-ray properties of each tar- 
get. We first give the number of Chandra obser- 
vations with >0 net source counts. We then show 
the line of sight absorption and photon index used 
for the model used to convert from intensity to lu- 
minosity for all HRC observations, and for ACIS 
observations with <200 net source photons. We 
finally give the best fit 0.3-10 keV unabsorbed lu- 
minosity and x^/dof for the best fit line of constant 
intensity. The lightcurves of the variable sources 
are presented in Fig. [TJ the complete version of 
this figure is available in the electronic version. 

3.1. Variability analysis 

We found significant X-ray variability in 28 of 
the 34 GCs and GC candidates; the remaining 
6 sources had 0.3-10 keV luminosities <2xl0'^^ 
erg s~^. The two AGN in our survey, B042D and 
B044D were rather faint (< 2 x 10^^ erg s~^), but 
both exhibited variability at at > 3a significance. 

In Fig. [5] we present a histogram showing the 
number of pairs of observation with a given sepa- 
ration T for XB144, the target with the most ob- 
servations. The histogram is binned logarithmi- 
cally, with a width of 0.2. Each bin averages tens 
or hundreds of observation pairs. 

We note that SF(r) becomes imaginary if the 
observed variability is smaller than the noise. Such 
behavior manifests as SF(r) = 0. Both AGN have 
mostly imaginary SFs; however, both show signifi- 
cant variability over times-scales of a few thousand 
days; this variability is consistent with t he ens em- 
ble AGN SF created bv lVagnetti et al.l (|201l[ ). 

The SFs of the brightest X-ray binaries (0.3-10 
keV luminosity >5xl0'^^ erg s~^) are also gener- 
ally consistent with the ensemble AGN SF, or even 
less variable, despite the relatively high signal to 
noise. This result is in keeping with the observed 



behavior of Galactic X-ray binaries; iMuno et al 



(|2002[ ) showed that the brightest X-ray binaries 
(i.e. Z-sources) vary by a factor of a few, while 



the fainter binaries (atoll sources) can vary by 1 
or 2 orders of magnitude. However, these bright 
X-ray sources have 2-10 keV fluxes greater than 
nearly all AGN; the 2-10 keV log N-log S plot for 
AGN calculated bv lMoretti et alj (|2003l) leads us 
to expect < 1 AGN per square degree with compa- 
rable fluxes. 

Our most encouraging result is that most of the 
fainter X-ray sources have SFs showing substan- 
tially more variability than the ensemble AGN SF 
over most time-scales. This is an important dis- 
tinction to make, since X-ray binary and AGN 
emission spectra are often very similar. 

Unlike the other bright X-ray sources, the SF 
for XB158 shows considerably more variability 
than the ensemble AGN SF over all time-scales. 
XB158 is a high inclination X-ray binary that 
exhibit s periodic intensity dips on a ^10,000 s 
period dTrudolvubov et al.ll2 002') and also disc pre- 
cession ( Barnard et al.l I2006D : it is further dis- 
cussed in Section [331 



3.2. Eight black hole candidates 

We have previously identified six black hole 
candidates (BHCs) in the observed region, based 
on their exhibition of characteristic low state emis- 
sion spectra at luminositie s that are too high 
for neutron star syst ei ns (IBarnard et al.l l2008t 
Barnard fc Kolbl 120091: iBarnard et al.l l2011aMbh . 
Five of these systems are associated with GCs, 
and four of those are persistently bright, which 
is unexpected given that much of the theoretical 
work to date models thes e systems as transients 
(see IBarnard et al.ll2011al for an overview of the 
literature on the topic of black holes in globular 
clusters). 

Low state emission spectra are characterized by 
a dominant power law com ponent with F <2, con- 
tribu ting >90% of the flux (jMcClintock fc RemillardI 
20061 ). and are seen in all low mass X-ray bina- 
ries, whether they conta in black holes or neutron 
stars (jvan der Klislll994l ). However, it has become 
apparent that the low sta te is confined to lumi- 
nosities 5,10% E ddington ([Gladstone et al.ll2007 : 
Tang et al.ll201ir) . Hence X-ray binaries that ex- 
hibit low state spectra at luminosities significantly 
higher than 3x10'^^ erg s~^ are likely to contain 
black hole accretors. 

However, a good fit with F <2 is not suffi- 



cient by itself to identify a BHC; we must first re- 
ject competing emission models. In particular, we 
must reject a disk blackbody-dominated spectrum 
(representing the thermally dominated black hole 
state), a blackbody -I- power law model that rep- 
resents high luminosity emission in neutron star 
XBs, and a disk blackbody -I- steep power law (F 
>2.4) model to represent the BH steep power law 
state. 

The Terzan-5 X-ray transient IGR J17480-2446 
exhibited an unusually impressive range of prop- 
erties: Z-source and atoll behavior, as well as 
therm onuclear X-ray bur s ts an d 11 Hz pulsa- 



tions; IChakraborty et al.l (|2011l ) conducted de- 



tailed spectral analysis of the outburst, using data 
from ^^40 ^daily RXTE observations. They found 
that the spectra were all well described by a sim- 
ple emission model (blackbody -I- power law -|- 
Gaussian emission line). The blackbody compo- 
nent contributed ~30-50% of the 3-15 keV flux, 
with kT = 1.4-2.1 keV for the most part and 3.7 
keV for the flrst observation; F ~2. 1-3.0; the un- 
absorbed 3-15 keV flux was 2.58-17.94x10"^ 

-2 t,-l 



erg 



cm ^ s ^ ( Chakrabortv et al.ll201l| ). 

We estimated the unabsorbed luminosity and 
blackbody contribution to each spectra in the 0.3- 
10 keV band for IGR J17480-2446, from the best 
fit parameters, and assuming a distance of 5.5 kpc. 
We assumed that the power law component ex- 
tended down to 0.3 keV; alternative models where 
the seed photon energy for the Comptonized com- 
ponent matched the inner disk temperature were 
rejected by the XMM-Newton spectrum of our 
M31 black hole candidate RX J0042.3+4115 i n 
the steep power law state (JBarnard et al.ll2011bl ). 
We found that the blackbody contributed 7-13% 
of the unabsorbed 0.3-10 keV luminosity of 1.8- 
33x10^'^ erg s'^ IGR J17480-2446 appears not 
to have been in the low state at any time during 
the outburst. 



We note that iLin et al.l (|2007l ) devised an al- 
ternative model to describe two other neutron 
star X-ray transients, Aql X-1, and 4U 1608— 52; 
traditional emission models (thermal -I- inverse- 
Comptonized emission) all yielded good fits for 
spectra with >10^ net source counts, but the au- 
thors raised objections to various aspects of the 
spectral evolution of the transient outbursts. They 
preferred a disk blackbody -I- blackbody emission 
model at high luminosities, with a Comptonized 



component necessary at lower luminosities. In 
this scenario, the neutron star transient evolution 
closely follows black hole transient evolution, with 
an extra blackbody component from either the 
neutron star itself or the boundary layer between 
the disc and neutron star surface. 

However, their model contradicts a large body 
of work based on observations of Galactic neu- 
tron star XBs. Plenty of evidence exists for an 
extended X-ray emission region in the high incli- 
nation, "dipping" X-ray binaries (radius ~20,000- 
700,000 obtained f rom ingress of photo-electric 
absor ption dips; see lChurch fc Balucihska-Church 
2004L and references within). Furthermore, broad- 
ened emission lines in Chandra observations of 
Cygnus X-2 su ggest a hot, dense corona of up to 



100,000 km (|Schulz et al.ll2009l ). Nevertheless, 



we applied the double-thermal emission model to 
our BHC spectra for completeness; the Comp- 
tonized component should not be requir ed at such 
high X-ray luminosities according to 
()2007l) . 



Lin et al. 



We examined the lightcurves and spectral histo- 
ries of all our targets, searching for further BHCs. 
We found a total of 8 GC LMXBs that apparently 
exhibited high luminosity low states; we refer to 
the X-ray source associated with GC labeled in 
the RBC v4 as Bnnn as XBnnn. In addition to 
the known BHCs XB082, XB144, XB153, XB163 
and X B185 (.Barnard fc Kolbll2009HBarnard et al 
2011al ). we have identified XB086, XB135 and 



XB148 as new BHCs. We present the lightcurves 
and spectral histories (F vs. time) for the 8 BHCs 
in FiglH the complete figure is available in the elec- 
tronic edition. The best fit blackbody -1- power 
law emission models for XB082, XB144, XB153, 
XB163, and XB185 are presented in Barnard et 
al. (2009, 2011); they clearly demonstrate the low 
state nature of the emission. 

In the following spectral analysis, uncertainties 
in best fit parameters and luminosities for XB086, 
XB135, and XB148 are quoted at the 90% confi- 
dence level, corresponding to ^1.6ct. 



3.2.1. XB086 



4.4 at 0.3-10 keV lumi- 
'^ . Furthermore, it ex- 



XB086 exhibited F 
nosities ^6-10x10^^ erg s~ 
hibited spectral variability ^^1600 days into the 
observations; F changed from 2.0 to 1.5 accom- 



panied by a luminosity decrease from 1.1x10 to 
7.5x10'^^ erg s~^. This variation may indicate a 
state change from e.g. the steep power law state to 
the low state. We note that XB086 apparently ex- 
hibited two different spectral states at very similar 
luminosities; hysteresis in the state transitions of 
black hole binaries where the transition from high 
state to low state is at a lower luminosity than 
the transition f rom low sta.te to high state is well 
known (see e.g. lMaccaronell2003[) . 

None of the ACIS spectra were of sufficient 
quality to discriminate between spectral states, 
so we obtained a spectrum for XB086 from the 
2002 January 6 XMM-Newton observation (ObsID 
0112570101). This spectrum rejected a disk black- 
body model (kTin = 1.7 keV xVdof =284/155), 
but a power law model provided a good fit (F = 
1.48±0.05, xVdof = 155/154). 

Adding a blackbody to the power law resulted 
in a slightly improved fit (kT = 0.9±0.3 keV, F 
= 1.47±0.18, xVdof =141/152); the blackbody 
component contributed ^--^12% of the 0.3-10 keV 
luminosity {^dtxlQ^'^ erg s~^). The Terzan 5 NS 
transient exhibit ed kT ~2.0 and F ~2.3 at sim- 
ilar luminosities ( Chakrabortv et al.l 12011 ): these 



parameters differ by >5cr. Fig. [5] shows the best 
fit two component model to the unfolded spec- 
trum, multiplied by channel energy; it is represen- 
tative of the BH low state, and inconsistent with a 
high luminosity NS XB. This spectrum containted 
■~9000 net source counts. 

Fitting a blackbody + disk blackbody model to 
B086 resulted in a 1.52J:o;^3 keV blackbody and an 
inner disk temperature of 0.66±0.08; the black- 
body component contributed 64% of the unab- 
sorbed flux. The inner disk temperature obtained 
for XB086 is > 5(t below th e mini mum tempera- 
tures obs e rved bv iLin et al.l (120071 ): furthermore, 
Lin et al.l (|2007l ) never observed blackbody contri- 



butions >50% in any of their transient spectra. 
Therefore the B086 spectr um is quite unli ke any 
of the spectra observed by iLin et al.l (|2007l) . 

We conclude that XB086 was in a low state; 
A^H — 1.06±0. 15x10^^ atom cm~^, and the unab- 
sorbed 0.3-10 keV luminosity was 6.4±0.2xl0^'^ 
erg s^^ (90% confidence limits). This result sup- 
ports our BHC classification from the ACIS obser- 
vations. 



3.2.2. XB135 

XB135 exhibited F '^1.6 at incredibly high lu- 
minosities: '-^4-6x10"^® erg s~^ during the Chan- 
dra observaations. The highest quality spectrum 
(ObsID 14198) had a ^^40 ks exposure time, and 
a count rate ~0.2 count s~^. This high intensity 
would result in significant pile-up if XB135 were 
on-axis; however, it is near the edge of the field of 
view, and the photons are distributed over a large 
number of pixels. Indeed, no pixel accumulated 
more than 19 photons over 40 ks, hence pile-up 
was not an issue. 

The ObsID 14198 0.3-7.0 keV XB135 emis- 
sion spectrum was well described by an absorbed 
power law (T = 1.72±0.08, xVdof = 205/221), a 
disk blackbody (kTin = 1.63±0.08 keV, xVdof = 
197/221, or a blackbody 4- power law model (kT 
= 0.7±0.2 keV, r = 1.5±0.4, xVdof = 187/219). 

We have previously published a nalysis of an 



XMM-Newton observation of XB135 ( Barnard et al. 



20081 ). It was well described by a power law emis- 
sion model (r = 1.56±0.03, xVdof = 467/435, but 
the fit was significantly improved by adding a ther- 
mal component (kT = 0.8±0.2, T = 1.54±0.14, 
X^/dof — 413/433); the thermal component con- 
tributed 11±5% of the flux. A disk blackbody- 
dominated model was rejected. We present the 
best fit blackbody 4- power law emission model to 
the pn spectrum in Fig. [51 this spectrum contains 
^15000 net source counts. 

The Terzan 5 transient exhibited kT 
keV and F 



2.76_g'Qg at its peak flux (|Chakrabortv et al 



1 fi4+°05 



20111 ): these are extremely different from the ob- 
served parameters for B135. Fitting a blackbody 
-I- disk blackbody model resulted in a l.Slg'g keV 
blackbody contributing 49% of the 0.3-10 keV un- 
absorbed luminosity (3.7x10'^® erg s~^); the inner 
disk temperature was 1.0±0.2 keV. Such a low 
disk temperature is consistent only with the low- 
est lumin o sities of the NS transients measured by 
Lin et al.l (J2Q07I ). and the blackbody contribution 
is too large. We conclude that this B135 spectrum 
is unlike either model for NS high states, and most 
liklely represents a low state BH XB. 

The inferred mass for XB135 is very large, at 
least 33 Mq from the XMM-Newton observation, 
and possibly >50 if XB135 was also in the low 
state during the ACIS observations. Such a mas- 
sive black hole could have been formed from a 



metal poor high mass star, so that little mass 
was lost in stellar winds during the aging pro - 
cess (|Belczvnski et al.ll2010l : iMapelh et al.lboiol l. 
or could have b een formed by stellar me rgers in 
the cluster core (JMiller fc Hamilton|[20Q2l ). 



3.2.3. XB148 

XB148 exhibited a hard emission spectrum (F 
-1.5) up to a 0.3-10 keV luminosity of 8x10^^ 
erg s^^. Unfortunately, none of the ACIS spec- 
tra suggesting high luminosity low states were of 
sufficient quality to discriminate between emission 
models. 

A 2008 January 27 XMM-Newton observation 
of XB148 (ObsID 0505720501) with -2000 net 
source counts yielded a 0.3-10 keV spectrum that 
was well described by a power law (F=1.80±0.15, 
xVdof = 96/115) or a disk blackbody (kTin = 
1.14±0.15, xVdof = 125/116); these resuhs are 
consistent with a black hole low state or high state. 
Fitting a blackbody -t- power law model resulted 
in kr = 0.4 keV, unexpectedly low for a neutron 
star system; we present this best fit to the unfolded 
spectrum times channel energy in Fig. [71 

The Terzan 5 transient exhibited kT ~ 1.4- 
3.7 (?), for 0.3-10 keV luminosities 2.5-33.7x10^7 
erg s~^, >5cr higher than for XB148; similarly, F 
for XB148 is >4cr lowe r than the range ovbserved 
by IChakrabortv et all ( 2011 ). Fitting a black- 
body -I- disk blackbody model yielded kT = 3.5 
keV, with unconstrained uncertainties, contribut- 
ing 42% of the unrealistic 0.3-10 keV luminosity 
(1.3x10'^^ erg s~^); the inner disk temperature is 
0.80ib0.1 0, 2>a below the lowest temperature mea- 
sured bv lLin et al.] (|2007f ). 

These results are supportive of the classification 
of XB148 as a black hole candidate. Assuming a 
power law model yields A^h = 1.5±0.4xl0^^ atom 
cm~^, giving an unabsorbed 0.3-10 keV luminos- 
ity of 4. 7±0. 5x10^7 erg s~^; hence, we cannot rule 
out a particularly massive (~2.8 M©) neutron star 
accretor in its low state. 

3.2.4. Properties of the host clusters 

We now consider the properties of the GCs 
hosting our BHCs. In particular, we examine the 
masses, ages, and metallicities of these clusters; 



these values were obtained from Caldwell et al. 



(|2011l) . who also ranked the 379 M31 GCs in their 



study by mass and by metallicity. We list these 
properties in Table |3l 

In I Barnard et al.l ( 2011al ) we concluded that 
particularly massive or metal rich GCs could con- 
tain bright X-ray sources (and hence BHCs) as 
well as those GCs that are both massive and metal 
rich. The GCs containing our three new BHCs 
are more massive than >80% of the GCs sur- 
veyed bv I Caldwell et al.l (|2Qlll ): however, two of 
them (B086 and B135) are in the lowest quar- 
tile for metaUicity, with [Fe/H] = -1.82. We 
note that the low metallicity for B135 is in keep- 
ing with the high infe rred BH mass for XB135; 



■■ f';-^ 

Belczynski et al.l (|2010l ) have shown that while BH 



masses are limited to ~15 M0 for solar metalic- 
ities, they could theoretically reach ^86 A/q for 
metallicities ^1% Solar. 

3.3. The dipping, precessing X-ray binary 
XB158, and related systems 

XB158 (a.k.a. Bo 158) is a high inclination 
XB that exper iences periodic intensity d ips on the 
orbital period (JTrudolvubov et al.ll2002l ) ; the orig- 
inal study used lightc urves folded on the 10017 
s (2.8 hr) period, and iTrudolvubov et al.l (|2002[ ) 
concluded that the dip depth varied with lumi- 
nosity. However, detailed analysis of the original 
XMM-Newton observations, along with propri- 
etary observations, revealed that the dips ap- 
pear ed in some ob servations, but not in oth- 
ers ( Barnard et al.l |2006). This behaviour oc- 
curs in systems with low mass ratios (short pe- 
riods) because the outer edge of the accretion disc 
reaches the 3:1 resonance, causing additional tidal 
torques from the secondary that lea d to elonga- 
tion and precession of the disc (see e.g. Osaki"l989|; 
Whitehurst fc Kindll99lHOgilvie fc Dub us 2001). 

XB158 is located at a high off-axis angle from 
the focus of the observations, and was only ob- 
served in 19 ACIS observations; it was observed in 
all of the HRC observations thanks to the larger 
field of view. The 0.3-10 keV luminosity was ex- 
tremely variable, ranging over ~4-40xlO'^^ erg 
s"^, with the best fit line of constant intensity 
yielding xVdof = 3030/63 for ACIS and HRC ob- 
servations. 

We note that the high quality XMM-Newton 
spectra of XB158 require a two component emis- 
sion model consisting of e.g. a black body -I- power 



law ( Barnard et al.l 120061) . however, the Chandra 
spectra were well described by a power law with 
best fit constant T ^ 0.59±0.17. 

We propose that the observed fluctuation be- 
tween high and low luminosities may be caused 
by varying accretion rates over the disc preces- 
sion cycle. Using proprietary three dimensional 
smoothed particle hydrodynamics code, we esti- 
mated the disc precessi on period to be ~29 o rbital 
periods, or ~81 hours ([Barnard et al.1120061) . The 
cadence of these observations is not well suited for 
testing such a period, however the similar minima 
and maxima are suggestive of some sort of super- 
orbital cycle. 

We note that the HRC luminosities for XB158 
may be misleading, since we are approximating 
a two-component spectrum with a single power 
law, and the HRC responses to each component 
may be different from the ACIS responses. How- 
ever, it is encouraging to see that ACIS and HRC 
observations made at similar times produced con- 
sistent luminosities. Even if the normalization of 
the HRC lightcurve is incorrect, the variability re- 
flects the intensity variation of the source, suggest- 
ing that the luminosity of XB158 really is rapidly, 
and systematically variable. 

We may expect GC XBs to have systematically 
shorter periods than XBs outside GCs because the 
probability of collisions with other stars in the GC 
is high, and these interactions shorten the orbital 
periods. Hence, even binaries with mass ratios 
>0.3 may have sufhciently short periods to expe- 
rience the tidal resonances that produce precessing 
discs. 

The 0.3-10 keV lightcurve of XB146 looks sim- 
ilar to that of XB158, in that it also fluctuates be- 
tween fairly consistent minima and maxima; how- 
ever, it is observed in 131 out of 134 observa- 
tions, and appears to be more structured. The 
best fit line of constant intensity yields x^/dof = 
3614/130, as the 0.3-10 keV luminosity fluctuates 
between ~ 10"^^ and ^SxlO'^^ ergs~^, with the oc- 
casional outlier. We therefore suggest that XB146, 
and other X-ray sources exhibiting similar behav- 
ior are XBs with precessing accretion discs. Un- 
fortunately, this is not easily verified unless they 
are viewed close to edge on. 



3.4. X-ray transients 

Four X-ray sources in our sample are consid- 
ered transient, in that they vary in 0.3-10 keV 
luminosity by a factor >100 between observations 
where they are detected, and are fainter than the 
detection limit in other observations. Two types 
of XB behaviour leads to transient X-ray sources. 
LMXBs may exhibit transient behaviour due to 
instabilities in the disc; the discs accumulate mat- 
ter during a long quiescent phase, and burn a large 
portion of the accreted material during outbursts. 
Long period HMXBs with eccentric orbits may 
also be X-ray transients, if accretion occurs only 
in a limited phase range near periastron. 

Three of the transients are associated with a 
GC or candidate and are expected to be LMXBs 
due to the ages of the GC populations; this in- 
cludes XB163, which we have identified as a black 
hole candidate, as previously discussed. 

XB128 varied by a factor >500 over 47 ob- 
servations, peaking at a 0.3-10 keV luminosity 
~5xl03^ erg 8"^ The best spectrum for XB128 
was from Chandra observation 4682 (2004 May 
23), yielding 157 net source counts; it was well 
described by a power law with r=1.5lQ;7, ^ ^^^^ 
blackbody with kT = 1.6to;4 keV, or a two compo- 
nent model consisting of a ~1 keV blackbody and 
a power law with T ~2; this spectrum is consistent 
with a neutron star or black hole primary. 

XPB-in7 varied by a factor >100 over 103 ob- 
servations. The peak luminosity of XPB-in7 was 
in Chandra observation 8183 (2007 January 14), 
with 227 net counts; however, we were unable 
to obtain an acceptable fit from any of our spec- 
tral models. Since the peak luminosity was only 
1.7x10'^^ erg s~^, we would be unable to deter- 
mine the nature of the accretor in any case. We 
note that PB-in7 is only a GC candidate, and fur- 
ther insight into the nature of the host may help 
interpret the nature of XPB-in7. 

XSK059A is coincident with likely compact GC, 
in a region rich in HII (Caldwell, private com- 
munication). As such, it may be a HMXB as- 
sociated with the HII instead of a LMXB asso- 
ciated with the GC. XSK059A varied by a fac- 
tor ~300 between Chandra observations, exhibit- 
ing long quiescent intervals punctuated with many 
outbursts. The lightcurve is consistent with 8 out- 
bursts >10'^^ erg s~^ separated by 1-6 cycles of 



~120 days, with the outbursts lasting up to ~90 
days; we indicate these bursts in Fig. 1 with down- 
ward arrows. However, these outbursts are not 
present in every cycle. Such behaviour may be ex- 
plained by a HMXB with variable mass loss rate 
from the secondary. 

4. Discussion and conclusions 

We have analyzed the long term lightcurves of 
globular cluster X-ray sources using 134 Chandra 
observations taken over '^12 years; these X-ray 
sources are likely LMXBs. Our sample was drawn 
from the Revised Bologna Catalogue V. 4, and in- 
cluded 30 confirmed old globular clusters, 4 candi- 
date globular clusters, 2 background galaxies, and 
1 star. Fitting each lightcurve with a line of con- 
stant luminosity revealed that 30 out 37 targets 
were variable; the other 7 targets had 0.3-10 keV 
luminosities <2 x 10'^^ erg s~^ . We infer from these 
results that all XBs in our field with 0.3-10 keV 
luminosities ^ 2 x 10'^^ erg s~^ should be variable 
also. 

We created structure functions from the 0.5- 
4.5 keV fiuxes of our variable sources, for com- 
parison with the ens emble AGN SF created by 
Vagnetti et al.l (|201ir) . The lower luminosity XBs 
exhibited SFs that generally showed substantially 
more variability than the AGN SF over most time- 
scales. The higher luminosity XBs had SFs gener- 
ally indicating similar, or lesser, variability than 
the ensemble AGN SF; however, the 2-10 keV 
fluxes for these systems were higher than almost 
all AGN. This discrimination is important, since 
the emission spectra of XBs and AGN are often 
similar. 

We have identified a total of 9 black hole candi- 
dates in M31 GCs, of which 8 are covered by this 
survey; 7 BHCs were identified in these data (4 re- 
ported in Barnard et al., 2011, and 3 new to this 
work), using our well established method of seek- 
ing low state emission spectra at luminosities that 
conspicuously exceed the limit for neutron star bi- 
naries. We have identified apparent super-orbital 
variability in XB158, which is a high inclination 
system where the disk is known to precess. The 
apparently regulated behaviour of XB146 may also 
indicate a short period system with a precessing 
disc. 

We will be creating SFs for the remaining ^^400 



X-ray sources in our field, and expect to identify 
>100 new X-ray binaries. 
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Fig. 1. — Long-term 0.3-10 keV luminosity lightcurves of each variable X-ray source, consisting of up to 89 

ACIS and 45 HRC observations. ACIS and HRC observations are indicated by red circles and black crosses 

respectively. For each lightcurve we provide the best fit line of constant intensity. At the top left of each 

plot we list the best fit constant intensity in units of 10^^ erg s~^ and the corresponding x^/dof. Lightcurves 

for the remaining GCs are presented in the electronic edition. 
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Fig. 1. — continued. The outbursts exhibited by SK059A are indicated by downward arrows. 
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Fig. 1. — continued 
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Fig. 2. — Histogram showing the number of observations with separation t for B144, which is the GC that 
was observed most often. 
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Fig. 3. — Struct ure functions for each of the variable X-ray sources, created from 0.5-4.5 keV fluxes for 
comparison with IVagnetti et al.l (J2011I ): the dashed line represents their ensemble AGN structure function. 
The dotted line shows ct^ for each bin for reference, since the noise component has already been subtracted 
from SF(t). We provide the name of the X-ray source at the top of each panel, along with the best fit 
constant luminosity. The y axis is scaled to the SF,^ and high values for the noise may be excluded. The 
complete figure is available online. 
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Fig. 3. — continued 
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Fig. 3. — continued 
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Fig. 3. — continued 
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Fig. 3. — continued 
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Fig. 4. — Long term 0.3-10 keV lightcurves of black hole candidates, along with the best fit photon index, 
r, for each ACIS observation; for observations with <200 net source counts T is fixed to the mean F for 
bright observations. Luminosity uncertainties are quoted at the la level, while F uncertainties are at the 
90% confidence level. The complete figure is provided in the electronic edition. 
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Fig. 4. — continued 



24 



Unfolded Spectrum 

XB086 
— I 1 — I — I — I — I — I 1 




I 
> 

CD 

7 

I 



o 7 

C/3 O t" 

o 
o 

> 

CD 



I 

o 



T 1 1 1 1 r 



HfPH 



0.5 




± 



.1 L. 



1 2 

Energy (keV) 



Fig. 5. — Unfolded XMM-Newton pn spectrum for XB086, assuming the best fit absorbed blackbody + 
power law model; the flux of each channel is multiplied by the energy. The power law component has similar 
r to the simple power law fit. 
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Fig. 6. — Unfolded XMM-Newton pn spectrum for XB135, assuming the best fit absorbed blackbody + 
power law model; the flux of each channel is multiplied by the energy. The power law component has similar 
r to the simple power law fit. 
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Fig. 7. — Unfolded XMM-Newton pn spectrum for XB148, assuming the best fit absorbed blackbody + 
power law model; the flux of each channel is multiplied by the energy. The power law component has similar 
r to the simple power law fit. 
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Table 1: Optical identifications for X-ray sources associated with members of the RBC V4. For each X-ray 
source, we provide th e X-ray position, along with the name of the object from the RBC. We then give the 
I Caldwell et al.l ( 2011 ) classification of the object, the distance between the X-ray source and its counterpart 



(Dobj), and the uncertainties in the X-ray position ( ctra and ctdcc)- Finally we give the deviation of the 
X-ray position from the center of the optical counterpart (Dist/cr). X-ray sources with no reliable centroid 
solution are indicated by ". . ." . 



Src 


X-ray Position 


Obj 


Class 


^Obj/" 


"■ra/" 


crDcc/" 


Dist/cr 


1 


00:41:37.598 -1-40:59:20.33 


SK081B 


Star 


2.6 








2 


00:42:6.011 +41:2:48.12 


B042D 


Gal 


0.9 


0.63 


0.55 


1.08 


3 


00:42:7.070 +41:0:17.55 


B044D 


Gal 


1.15 


0.64 


0.67 


1.24 


4 


00:42:12.198 +41:17:58.78 


B078 


OldGC 


0.54 


0.28 


0.27 


1.41 


5 


00:42:15.847+41:1:14.98 


B082 


OldGC 


0.68 


0.41 


0.37 


1.23 


6 


00:42:18.670 +41:14:02.06 


B086 


OldGC 


0.2 


0.27 


0.27 


0.52 


7 


00:42:19.769 +41:21:53.99 


r3-74 


PosGC 


0.55 


0.33 


0.37 


1.12 


8 


00:42:21.588 +41:14:19.92 


SK050A 


Pos GC 


0.38 


0.28 


0.27 


0.99 


9 


00:42:25.034 +40:57:19.03 


B094 


OldGC 


1.43 


3.0^^ 


2.8" 


0.4 


10 


00:42:26.074 +41:19:15.02 


B096 


OldGC 


0.23 


0.27 


0.27 


0.60 


11 


00:42:27.392 +40:59:36.56 


B098 


OldGC 


0.69 


3.1" 


2.5" 


0.2 


12 


00:42:31.271 +41:19:38.95 


B107 


OldGC 


0.24 


0.27 


0.27 


0.63 


13 


00:42:33.161 +41:03:28.09 


BllO 


OldGC 


0.76 


1.6^* 


1.2" 


0.4 


14 


00:42:34.311 +40:57:08.86 


B117 


OldGC 


0.66 


2.2° 


1.4" 


0.3 


15 


00:42:40.664 +41:10:33.67 


B123 


OldGC 


0.38 


0.28 


0.28 


0.96 


16 


00:42:41.451 +41:15:23.93 


B124 


OldGC 


0.26 


0.27 


0.27 


0.68 


17 


00:42:46.101 +41:17:36.35 


PB-in7 


PosGC 


0.37 


0.39 


0.36 


0.70 


18 


00:42:47.814+41:11:13.94 


B128 


OldGC 


0.22 


0.27 


0.27 


0.58 


19 


00:42:50.777 +41:10:33.80 


BH18 


OldGC 


0.71 


0.27 


0.27 


1.86 


20 


00:42:52.131 +41:31:05.40 


B135 


OldGC 


0.56 


0.42 


0.42 


0.94 


21 


00:42:55.623 +41:18:35.10 


B138 


OldGC 


0.33 


0.61 


0.53 


0.41 


22 


00:42:59.665 +41:19:19.25 


B143 


OldGC 


0.16 


0.28 


0.27 


0.41 


23 


00:42:59.872 +41:16:05.76 


B144 


OldGC 


0.23 


0.27 


0.27 


0.60 


24 


00:43:1.346 +41:30:17.04 


B091D 


OldGC 


1.08 


0.59"^ 


0.61" 


1.3 


25 


00:43:2.938 +41:15:22.60 


B146 


OldGC 


0.21 


0.27 


0.27 


0.55 


26 


00:43:3.309 +41:21:21.71 


B147 


OldGC 


0.21 


0.27 


0.27 


0.55 


27 


00:43:3.868 +41:18:04.87 


B148 


OldGC 


0.16 


0.35 


0.39 


0.31 


28 


00:43:7.526 +41:20:19.48 


B150 


OldGC 


0.07 


0.44 


0.52 


0.10 


29 


00:43:9.853 +41:19:00.67 


SK059A 


HII 


0.08 


0.43 


0.46 


0.13 


30 


00:43:10.622 +41:14:51.37 


B153 


OldGC 


0.06 


0.60 


0.64 


0.07 


31 


00:43:14.388 +41:7:21.42 


B158 


OldGC 


0.28 


1.17 


1.35 


0.16 


32 


00:43:14.535 +41:25:12.88 


B159 


OldGC 


1.4 


2.0" 


1.7" 


0.5 


33 


00:43:15.465 +41:11:25.51 


B161 


OldGC 


0.44 


0.59 


0.53 


0.56 


34 


00:43:17.710 +41:27:45.53 


B163 


OldGC 


1.0 


0.31 


0.31 


2.26 


35 


00:43:36.798 +41:08:13.46 


B182 


OldGC 


1.9 


2.9" 


3.7" 


0.11 


36 


00:43:37.309 +41:14:43.11 


B185 


OldGC 


0.53 


0.76 


0.96 


0.43 


37 


00:43:56.315 +41:22:01.16 


B204 


OldGC 


2.09 


1.1" 


1.3" 


1.2 



Positional data obtained from a heavily binned image where each pixel represents 9x9 native pixels. 
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Table 2: For each object, we provide the number of Chandra observations that include the object (ACTS 
and HRC). We then give the absorption (-/Vh) and photon index (F) for the model used to convert from 
source intensity to luminosity, followed by the x^/dof for A^h and F. Finally we show the unabsorbed 0.3-10 
keV luminosity and x^/dof for the best fit line of constant intensity. Uncertainties are quoted at the 90% 
confidence level. 



Src 


A"acis 


A^HRC 


Ah / 1020 


F 


xVdof (A^h) 


xVdof (F) 


r37 
-'^con 


xVdof {Lil,) 


1 


2 


1 


7 


1.7 






0.06 ± 0.06 


1.9/2 


2 


16 


31 


7 


1.7 






0.10 ± 0.02 


83/46 


3 


12 


5 


7 


1.7 






0.05 ± 0.02 


62/16 


4 


78 


45 


44 ± 12 


1.83 ± 0.18 


0.44 / 4 


1.5 / 4 


1.48 ± 0.03 


2455/122 


5 


17 


41 


50 ± 10 


1.5 ± 0.14 


5/16 


5/16 


19.7 ± 0.4 


172/57 


6 


76 


45 


17 ±7 


1.43 ± 0.1 


19/ 15 


6/15 


5.79 ± 0.07 


1490/121 


7 


70 


19 


7 


1.7 






0.055 ± 0.007 


101/88 


8 


79 


39 


21 ± 11 


2.2 ± 0.3 






0.257 ±0.013 


766/117 


9 


5 


6 


7 


1.7 






0.12 ± 0.03 


24/10 


10 


83 


45 


42 ± 12 


1.9 ±0.19 


0.6/3 


4/3 


2.2 ± 0.04 


2695/127 


11 


12 


19 


7 


1.7 






0.05 ± 0.02 


33/30 


12 


84 


44 


27 ±9 


1.81 ± 0.18 


3/3 


0.09 / 3 


2.37 ± 0.04 


3142/127 


13 


19 


29 


7 


1.7 






0.15 ± 0.02 


74/47 


14 


11 


10 


7 


1.7 






0.21 ± 0.03 


134/20 


15 


77 


43 


7 


1.7 






0.135 ± 0.010 


312/119 


16 


85 


44 


7 


1.7 






0.36 ± 0.14 


486/128 


17 


75 


28 


7 


1.7 






0.082 ± 0.007 


694/102 


18 


32 


15 


7 


1.7 






0.019 ± 0.005 


612/46 


19 


73 


26 


7 


1.7 






0.069 ± 0.007 


179/97 


20 


14 


45 


33 ±2 


1.56 ± 0.04 


66/ 13 


50 / 13 


48.9 ± 0.4 


348/58 


21 


85 


43 


7 


1.7 






0.281 ± 0.011 


759/127 


22 


83 


45 


16 ± 4 


1.85 ± 0.12 


1.4 / 11 


3/11 


5.03 ± 0.05 


483/127 


23 


86 


45 


9±4 


1.36 ± 0.07 


6/ 13 


7/13 


5.63 ± 0.06 


1310/130 


24 


9 





16 ± 4 


0.74 ± 0.12 


1.8/ 1 


1.1/1 


1.09 ± 0.09 


379/8 


25 


86 


45 


16 ± 5 


2.04 ± 0.13 


8/13 


4/13 


3.17 ± 0.04 


3614/130 


26 


81 


45 


30 ± 30 


2 ±0.4 






1.77 ± 0.03 


334/125 


27 


85 


45 


17 ± 4 


1.76 ± 0.12 


3/15 


12 / 15 


4.05 ± 0.05 


2568/129 


28 


82 


42 


7 


1.7 






0.371 ± 0.014 


474/123 


29 


77 


39 


27 ± 16 


1.9 ±0.3 


0.9/3 


3/3 


0.190 ± 0.010 


2361/115 


30 


83 


45 


12 ± 1.7 


1.51 ± 0.03 


24/ 79 


35 /79 


13.5 ± 0.12 


825/127 


31 


19 


45 


26 ± 15 


0.59 ± 0.17 


2.4/8 


4/8 


15.5 ± 0.3 


3030/63 


32 


23 


24 


7 


1.7 






0.069 ± 0.011 


48/46 


33 


81 


37 


7 


1.7 






0.166 ± 0.011 


252/117 


34 


15 





12 ± 12 


1.6 ± 0.3 


0.11 /2 


0.12/2 


0.11 ± 0.02 


428/15 


35 


9 


20 


7 


1.7 






0.049 ± 0.010 


26/28 


36 


30 


45 


11 ± 4 


1.51 ± 0.08 


10/25 


13/25 


6.81 ± 0.10 


196/74 


37 


31 


45 


7 


1.7 






0.84 ± 0.04 


284/75 
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Table 3: Ages, masse s and metallicities (expressed as [Fe/H]) of the GCs hosti ng our BHCs, taken f rom 
Caldwell et ail (|201l[ ). We also rank the GCs amongst the 379 CCs studied by ICaldwell et all (|201l[ ) by 
mass and metallicity. 



GC Log Age/ yr Log Mass/M© Mass Rank 



[Fc/H] [Fc/H] Rank 



B045 


10.2 


B082 


10.1 


B086 


10.2 


B135 


10.2 


B144 


10.1 


B148 


10.2 


B153 


10.1 


B163 


10.1 


B185 


10.1 



5.94 


69 (>82%) 


-0.86±0.11 


137 (>64%) 


6.72 


4 (>99%) 


-0.74±0.11 


102 (>73%) 


6.13 


33 (>91%) 


-1.82±0.16 


294 (>22%) 


5.97 


64 (>83%) 


-1.82±0.13 


296 (>22%) 


5.58 


139 (>63%) 


+0.08±0.17 


10 (>97%) 


6.04 


50 (>87%) 


-1.02±0.11 


181 (>52%) 


5.81 


87 (>77%) 


-0.28±0.13 


36 (>91%) 


6.26 


20 (>95%) 


-0.13±0.11 


21(>94%) 


6.03 


53 (>86%) 


-0.61±0.11 


75 (>80%) 
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